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Abstract. We show that, within some modified gravity theories, such as the 
Palatini models, the non-linear nature of the field equations implies that the 
usual naive averaging procedure (replacing the microscopic energy-momentum by 
its cosmological average) could be invalid. As a consequence, the relative motion 
of particles in Palatini theories is actually indistinguishable from that predicted 
by General Relativity. Moreover, there is no WEP violation. Our new and most 
important result is that the cosmology and astrophysics, or put more generally, 
the behaviours on macroscopic scales, predicted by these two theories are the 
same, and as a result the naturalness problems associated with the cosmological 
constant are not alleviated. Palatini gravity does however predict alterations to 
the internal structure of particles and the particle physics laws, e.g., corrections 
to the hydrogen energy levels. Measurements of which place strong constraints 
on the properties of viable Palatini gravities. 
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Extensions of General Relativity (GR) have always received a great deal of 
attention. Such theories are motivated by quantum gravity models and by the wish 
to find phenomenological alternatives to the standard paradigm of dark matter and 
dark energy [T]. 

Modifications to Einstein gravity generally result in non-linear (in energy 
momentum tensor T^i/) corrections to the field equations. The application of these 
equations equations to macroscopic {e.g. cosmological) scales involves an implicit 
coarse-graining over the microscopic structure of particles. However, when there are 
extra non-linear terms in the field equations, a priori, the validity of the usual coarse- 
graining procedure can no longer be taken for granted, as first proposed by [2J. Hence, 
as discussed in [2] , it may be important to take into account the microscopic structure 
of matter when applying the field equations to macroscopic scales. 

Unfortunately, to date, that has not been the common practice j3j. This is 
probably because, in GR, as in Newtonian gravity, the microscopic structure of matter 
is not particularly important on macroscopic scales. It is standard practice to replace 
the metric, g^,^, and the energy momentum tensor, T^^, with some average of them 
that coarse-grains over the microscopic structure. This simple procedure only works, 
however, because on microscopic scales the equations of GR are approximately linear. 

In this Letter we show that such an approach cannot simply be applied to modified 
gravity without a detailed first analysis of the energy-momentum microstructure. 
Indeed, naively averaging over the microscopic structure will generally lead one to 
make incorrect predictions, and inaccurate conclusions as to the validity of the theory 
[2]. Indeed, it is possible that a theory deviates significantly from GR at the level of 
the microscopic field equations to be indistinguishable from the latter when correctly 
coarse grained over macroscopic, e.g., cosmological, scales. 

We illustrate this point for a class of modified gravity theories in which the Ricci 
scalar, R, in the Einstein-Hilbert action is replaced by some function f{R,R^^R^^), 
i.e.: 



The field equations for this action can be derived according to two inequivalent 
variational approaches: metric and Palatini. In the former, R^^i, and R are taken 
to be constructed from the matter metric g^^, which couples to matter and governs 
the conservation of energy momentum tensor, and the field equations are found by 
minimizing the action with respect to variations in g^u- 

In the alternative. Palatini 

approach, R = R^^g^^ where R^i, is a function of some connection field which is, 
a priori, treated as being independent of g^j/. The field equations are then found by 
minimizing the action with respect to both V^^ and g^^. If f{R, R'^'^R^j^u) — R~ 2A 
{i.e. GR with a cosmological constant) then the two approaches result in the same 
field equations. Otherwise they are generally different. Note that the Palatini f{R) 
field equations are mathematically equivalent to an a; = — | Brans-Dicke theory with 
a potential, and thus represents more general modified gravity theories. 

Within the metric approach to the f{R) gravity theories, averaging over 
microscopic scales is generally no less straightforward than it is in GR; this is because 
in both cases all degrees of freedom are dynamical. These dynamics normally ensure 
that the field equations, for all degrees of freedom, are approximately linear (in 
energy momentum tensor T^i/) on small-scale structures [4]. In contrast, averaging 
in Palatini models is not so trivial as the new degree of freedom is non- dynamical, and 
so its field equation remains non-linear (in energy momentum tensor T)ii/) even on the 
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smallest scales. Consequently the cosniological behaviour of these theories can be very 
different from what has been suggested in the literature [21 [3 H] , which adopt the same 
averaging procedure as in the metric approach. Also, although GR plus a cosmological 
constant is a special case of Palatini f{R) (and also of metric f{RY) theories, the 
averaging problem for general Palatini theories dose not arise there because the GR 
field equations are linear (and algebraic) in curvature R and energy momentum tensor 
components, while in contrast for general Palatini theories the curvature depends 
nonlinear ly (but also algebraically) on energy momentum tensor. We emphasize again 
that it is this " nonlinearity" that makes the averaging of general Palatini theories not 
as trivial as that in GR. 

Palatini f{R^^'^R^^) theories are similar to the f{R) ones in many aspects, but 
their study is more complicated. In what follows we shall mainly focus on the latter 
and state the results of the former when appropriate, referring more details to [4]. To 
avoid confusion later, we now take i? — > 7?., where TZ = Rfi„{T)g'^'' . 

Varying the f{TZ) action with respect to gives that T^^ is the Levi-Civita 
connection of g^i/ = f'(JZ)gf^^. Varying the action with respect to g^i, gives: 

1 

kT'^^ =kT ^ /'($)$ -2/($), (2) 
where we have defined ^ = TZ, R = R^i^g^"" = $//'($), G^_,u{g) = - \9ij.vR and: 



G^(5) = 77i7:;:T['^^^-n$)^M, (i) 



vm = \ (/'($)$ - /($)) , T% = - ^^^p- . 

Also V^^V/^^ = 0, where v|f^g^p = defines V^F\ In f{R'"'R^,^) gravity the field 
equation is similar to Eq. ([1]), but the metrics and 5^1^ are disformal, and V 
depends algebraically on 7^^ rather than simply on T [4j. 

It is important to bear in mind that Eqs. ([T]) and ([2|) are microscopic field 
equations, and they are only definitely valid when all the microscopic structure 
in the distribution of energy and momentum is taken into account. Nonetheless, 
the cosmological and astrophysical behaviours of these theories have, to date, been 
studied by simply replacing 5^^, $ and Tf^^ in Eqs. (HI [J) by some coarse-grained 
averages of them [S1[31|S]. However, in Palatini f{TZ) theories, a deviation from GR 
requires that / depends nonlinearly on This nonlinearity introduces an averaging 
problem for the Palatini theories. As we shall show, this means that the standard 
averaging procedure is no longer valid and generally results in incorrect physical and 
mathematical predictions. Furthermore, when the microscopic structure of matter is 
taken into account, the late-time cosmology of essentially all f{TZ) Palatini theories 
is indistinguishable from that of standard GR with a cosmological constant. We 
also show that these these theories do not violate the weak equivalence principle 
(WEP). Palatini f{R) gravity is however a different theory from GR and does produce 
alterations to microscopic particle physics that could be detected. 

The microscopic structure of the space-time distribution of matter energy density, 
p, will not affect the macroscopic, or coarse-grained, behaviour of the theory if and 
only if the field equations of the theory are linear in p. Here is a simple example: 
consider a region of space with average density (p) and volume V, which contains N 
particles each with density pc and volume Vp. The space in between the particles is 
empty and so (p) V — NpcVp. Now consider some quantity Q{p). Inside the particles. 
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p = Pc and so Q — Qc = Q{pc)', outside p — and so Q ^ Qo = Q{0). The average 
value of Q (by volume) is: 



(Q) - Qo 1 



NVr. 



Qii 



Qc — Qo 



(p) 



(3) 



It is clear that, irrespective of how Q depends on p, (Q) depends on (p) linearly. Thus 
if Q depends nonlinearly on p, we have {Q{p)) ^ Q {{p))- In Palatini theories, V{^) 
generally exhibits highly non-linear dependence on p, and so {V{^{p))) ^ V {^{{p))) . 

It is also interesting to know when naively coarse-grained equations are valid 
to a good approximation over the length scales and density scales where we have 
observations. This only occurs in Palatini theories when the inherently nonlinear 
modifications to GR have sub- leading effect. In all these cases the theory will reduce 
to GR up to the order where the naively coarse-grained equations are valid. In this 
work we are concerned with the how Palatini theories deviate from GR, even when 
those deviations are sub-leading order. We must then always take account of those 
nonlinear terms in the equations which determine how Palatini theories deviate from 
GR. The crucial role played by these nonlinear terms means that the naive coarse- 
graining procedure will always fail to accurately describe the differences between 
Palatini theories and GR. 

To uncover how Palatini gravities behave on macroscopic scales, we consider how 
a set of microscopic particles evolve under such a modified gravity. Consider a single, 
spherical particle for which 7^^ ^ for R < Rp but vanishes otherwise. The metric 
for such a particle is [5]: 



1 



W{r) 



dr"^ +r'^dn'^, 



with 



W{r) = I - 



2GM{r) y($o) 



2GM{r) ^ K dr' r 
Jo 



r 



3/'2(<i>o) ' 
p AV{^) 



X 



, , , {p+Pi+Pa) 



(4) 
(5) 
(6) 



where 



AV^ = V{^) - /'2(ci>)l/($o)//'2(ci>o), 

- (P + Pi) u^'uPgp, + PI5^', + tt^,, 

u^^u'^g^^i, — —1 and tt^'j, is the anisotropic stress satisfying ir^^yu" = 0. In the rest 
frame of the particle 



= diag(0,pA, -pa/2, -pa/2), 



$0 = ^{T = 0). 



The presentation of this solution in Ref. [5] was, however, incomplete as it failed to 
note that the 69 component of Eq. JT]) results in the condition: 

d 



r — Pcff 

dr 



P + Pl ^PA 

f'^{^)rW{r) 



Y{r) 



PA 



where 



Y{r) = (^4^GPeffr' 



GM 



Aoffr^ 



(7) 
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Ks = Vi^o)/f'H%) and 

PeS^ipi+PA-^V/K.)/f'\ 

Eq. ([7]) implies that aX R — Rp one has Pcff — 0. 
Outside the particle, x = and 

W^(r) = 1 - GMp/r - AeftrVS 

where Mp = M{Rp) = const and so the metric is precisely that of a Schwarzschild 
de-Sitter spacetime with gravitational mass Mp and an effective cosmological constant 
Aeff. A similar Schwarzschild de-Sitter solution is found for the Palatini f{R^'^Rp.,j) 
case [4] , with V replaced by a more complicated function of T'^j^ . Thus we see that in 
Palatini models the external metric of a single particle is precisely what it would be 
in General Relativity with a cosmological constant. 

We now consider a spacetime containing many such particles. We define X(/)(t) 
to be the position of a particle / and w^^^ = i^^^ . Making a weak field approximation 
with respect to (/^i^ and using the above solution, we find that inside the particle 
labeled K and to C(e^), where e ^ max(|u(x)|, VTl — W{r))): 

g^.Ax^dx" = - {W{r) + 2x{r) - 2U{x)) At^ + (f + 2U{x)) dx^Ax^ 
Uix) ^ ^"^^^'^ 



^■^(if) ~ ^ ~ ''^A(i) is the active gravitational mass of each 

particle given by M{Rp). 

Let u^ji^ be the 4- velocity of the ii''^ particle which satisfies ~ 
To the order 0(e), 

Now T^i, is conserved with respect to i.e. vIF't^i/ = 0. Evaluating this equation 



at the centre of the particle we have u'^Wjf u'^ = 0, i.e 
dT2 ~ 2" 



d'xlj,) 1 



= -5'^3oo,,"''w" = -r{m,^■ (8) 



where r is the proper time along the worldline of the particle, and so dt/dr — f ^/2($) 
(because the metric — dr^ = + dx^) means 1 = -p- [(dt/dr)^ + {dx/dr)'^] w 

■^{dt/drf). 

The internal configuration of the particle can be static or non-static, and here for 
simplicity we assume it to be static and will comment on the non-static case below. 
This requires that all gradients in $ cancel with gradients in the pressure and also 
that T and hence <& are conserved along particle worldlines, i.e., u'^d^T = u^d^^ = 0. 
Equation ([5]) is then equivalent to: 

"(^)"^ = -^-"— 
The acceleration with respect to r depends both on the gravitational field, U^i, and 
/'($). The acceleration measured with respect to however, depends only on U^i. 



Indistinguishable Macroscopic Behaviour of Palatini Gravities and General Relativity^ 



The relative acceleration of two particles labeled 1 and 2 as measured by a third, 
labeled 3 say, with proper time r(3) is therefore: 

where dt/dT(^j-j — /'^/^($(3)). In a uniform gravitational field: U^i = const and so 

= 0(2) =^ A«l2 = 0. 

It follows that an observer {i.e. particle 3) sees any two other particles, 1 and 2, 
accelerate at the same rate in a uniform gravitational field. This is precisely what is 
required by the Weak Equivalence Principle. 

Since the internal configuration of the particles is static and the centres of two 
particles do not have relative acceleration, there will be no relative acceleration 
between the two whole particles. In Ref. [5] it was suggested that internal gradients 
in would lead to WEP violations. In fact, because those gradients all vanish outside 
the particle, they cannot affect the overall motion of the particle [8j (this is related 
to the fact that the extra force, which arises due to the modification to GR, has zero 
range; as we shall discuss in more detail below). In this case, hydrostatic equilibrium, 
Eq. ([7]), ensures that gradients in $ are cancelled by pressure gradients. 

The absence of WEP violation in the dynamics of particles ensures that the 
inertial and passive gravitational mass of particles are equal. Moreover, the inertial 
mass and the active gravitational mass of particles are also equal in these theories. 
Let us define 

'^T^cff. = G^(g)+Aeff5^. (10) 
By the modified Einstein equation, Eq. ([T]), and the contracted Bianchi identity: 

V^^^T'^offi/ — where '^'f^gi.p — 0. Outside an isolated particle T'^effi/ = and, as we 
have seen above, the metric g^^jy outside this quasi-static (i.e. Wj <C 1) isolated particle 
is Schwarzschild-de-Sitter. In Ref. [9], Tolman shows that, under these conditions, 
the inertial and active gravitational mass of an isolated (not necessarily spherically 
symmetric) system are equal (see also [8]). For our particles, the inertial, passive and 
active gravitational mass are therefore equal to nip where by Eq. ([5]): 



An r dr 



Kp + AV 



This equivalence holds for particles both in GR and in Palatini f{TZ, R^^^R^^) theories 
[4] . It follows that the motion of isolated particles (situated in a vacuum) in Palatini 
theories is exactly the same as it is in GR with A — Aoff. This has a very 
simple physical explanation. In all modified gravity theories one may think of the 
modifications as being due to some new, effective force. In Palatini theories, there is no 
extra dynamical degree of freedom, so this 'new' force is also non-dynamical, i.e., does 
not propagate and only acts at points. The closest analogue to this in particle physics 
would be Fermi's original proposal for a theory of the weak interaction. The effective 
new force is entirely local and depends on gradients in the T or T^. There is therefore 
no way for one isolated particle to influence the motion of another, and the force 
simply vanishes in a vacuum. The effective force will alter the internal configuration 
of the particle nonetheless, but this will not effect the motion of the particle as a whole. 
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Note that the above discussion reUes on the specific model for the matter distribution 
in a system. In most systems it is realistic to take the matter as to be distributed in 
localized clumps (what we refer to as particles) which only interact with each other 
gravitationally (or through extra forces due to the modification of gravity, which as we 
said above is of zero range and thus completely irrelevant for classical particles). For 
the extreme environments like a neutron star, however, the energy density is of order 
nuclear density, and there is essentially no space between different particles so that the 
above matter model breaks down: in this case we expect the naive averaging procedure 
to be applicable. For cosmology and most other astrophysical systems our averaging 
is nonetheless more realistic. We emphasize that the absence of WEP violation does 
not require the internal configuration of the particles is static (as we have assumed). 
Rather, it is a result of the fact that the extra force has zero range so that there is no 
composition-dependent extra force between different particles which do not overlap 
with each other. Additionally Palatini f{R) theories are equivalent to a generalized 
Brans-Dicke, and, in common with other Brans-Dicke theories, the extra force has no 
composition dependence. 

We now apply our results to a cosmological setting. At late times and on a 
microscopic level, most matter in the Universe is made up of small particles and a bit 
of radiation. In f{TZ) theories, T = for radiation and so the presence of radiation does 
not alter the relation between $ and T. Our analysis is therefore directly applicable 
to this setting. The Universe must therefore evolve precisely as it would in GR with 
a cosmological constant Aoff. This can be seen in an alternative manner: averaging 
Eq. (jlOp over volume V containing N particles with mass rrip and using Eq. ([3|) we 
find 

Using Eq. ([7]) we find that in the rest frame of the particle (^kTI^^^ = 0. More 

generally then I^kTI^^^ ^ 0{Kp'^^^^^j.Sv^6vj) where 5v^ is the relative particle 

velocity. Thus, when correctly coarse-grained over cosmological scales, T^^g,^ describes 
a coUisionless dust; we have assumed that the peculiar velocities, 6v, of the particles 
are small and dropped terms of 0{Sv'^). The cosmological evolution of the Universe 
in such theories is therefore precisely the same as it is in GR with a cosmological 
constant Aoff and a dust with energy density pJi^ttor- This argument also holds for 
most astrophysical systems. In particular, solar system tests are evaded and the 
Parametrized Post-Newtonian parameters are indistinguishable from those of GR. 

Many of the predictions made in the literature (e.g. 3 ) do not, therefore, follow 
from a Palatini f{TZ) model, but may still be correct for some other modified gravity 
theory. Any such theory would, however, likely be subject to additional constraints 
from local tests of gravity. 

Despite of their similar behaviours. Palatini f{Ti) gravities and GR are not 
equivalent. Even though motions of isolated bodies are the same in both theories, the 
internal structure and dynamics of them are generally different [3]. This is because, 
for the body to be stable under gravity one must require that Eq. ([7]) holds, and 
the appearance of AV in this equation clearly indicates an alteration to GR. As we 
mentioned above, the pressure gradients inside the particle must be modified with 
respect to those in GR in order to exactly cancel the gradients in <&, and so we 
expect the distributions of matter inside this particle, that is, the internal structures, 
according to these two theories to be different. 
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As another example, the electron energy levels in atoms are altered in Palatini 
fiTV) theories (see also [IT] for a recent work). Such alterations occur, because in these 
theories the effective electron mass depends on the electron density, and so is different 
for different energy levels. In the absence of any Palatini modification, we define the 
energy and number density of an electron with total angular momentum j , in energy 
level n, to be: E^j = —a^m^Snj and ni""''' respectively where a is the fine structure 
constant, and rrif, the electron mass. In Palatini f{TV) theories, the electron mass 
depends on $ and hence on the electron density which is different for different values 
of n and j. We assume that /(<&) « 6<f'(l + £($)) so that the density dependence of 
the electron mass is slight. The total energy of an electron (up to an overall constant) 
is then g]: 

where 

for some and to leading order in a^: 

A„, ^ - (e(<i>))„^. ^-J d=^xn("^)(x)e(|.„,(x)) (12) 

where $„j(x)/& = Kmgni""''' (x). Thus, in Palatini theories determinations of the 
electron mass from transitions between different energy levels using the standard 
formula for Eni would find different answers for each transition unless A„; = const; 
this possibility is however very strongly constrained. So far both e and rUe have only 
been defined up to an overall constant i.e. e —^ e + So, rrie rne{l + Sq) for some 
constant Jo *C 1 is allowed. We fix the definition of mg, and hence also e, so m^, is the 
effective electron mass for the ground state i.e. {s)^q = 0. Using measurements of the 
electron mass from the transitions 1S-2S and 2S-8D [12] we have 

< 8 X 10"^^. (13) 

Respectively for (10), (20) and (83) states the values of pe = nieUe near the expected 
electron radius are: 3 x 10~^gcm~^, 10~^gcm~'^ and 10~^gcm~^. If, for instance, 
e(<i>) w const + eo^/bHQ, where is the value of the cosmological constant today, 
Eq. p3p gives the very strong constraint: 

|eo|«|/"(<i>)i?oV/'('i>)l< 4x10-40 

In summary, much of our intuition about how the microscopic behaviour of 
gravity affects physics on large scales is based upon Einstein's general relativity. In 
this Letter we show that such an intuition cannot simply be generalized to modified 
gravity theories without a detailed analysis of the energy-momentum microstructure. 
Indeed, naively averaging over the microscopic structure will generally lead to incorrect 
predictions, and inaccurate conclusions as to the validity of the theory. In particular, 
the naive averaging procedure is invalid in Palatini theories. A correct averaging 
procedure shows that the cosmology of Palatini f{R) models is identical to that 
of GR and fine tuning problems associated with the cosmological constant are 
neither alleviated nor, it should be said, worsened. Furthermore, the relative motion 
of particles in Palatini theories is indistinguishable from that predicted by GR. 
Interestingly, although Palatini f{R) theories were designed to modify gravity on large 
scales, they actually modify physics on the smallest scales {e.g. the energy levels of 
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electrons) leaving the larger scales practically unaltered. In general, before considering 
any astrophysical consequences of a modified gravity theory, it is important then to 
check that it does not make unrealistic predictions for atomic physics. 

One may wonder whether similar problems arise in the metric f{R) gravity 
theories. In this paper, we have been concerned with averaging over very small scales 
(e.g. atomic scales). In the metric f{R) theories, the extra scalar degree of freedom 
is dynamical and so it can and does propagate. This means that its dependence on 
the distribution of matter is not so rigid as in Palatini theories. Over the very small 
scales of particles, although the energy density of matter might change rapidly (from 
nuclear density inside the particles to zero outside), the extra scalar degree of freedom 
$ is not required to change so abruptly. In metric f{R) theories, the field equation 
for $ has the form: + V'{^) oc (p — 3p), where p is the pressure of matter and p 
its energy density. Over very small scales, the kinetic term (V^^) dominates over the 
potential term {V') and so reduces to V^fE" cx (p — 3p). This means that over small 
scales the leading order deviation from GR, determined by $, obeys a linear second 
order differential equation. In Palatini theories, $ is related algebraically to p — 3p. 
If $ depends algebraically and linearly on p — 3p, we reduce to GR. It follows that in 
Palatini theories the leading order deviation from GR is necessarily determined by a 
non-linear algebraic relation to p — 3p. 

One may think of metric f{R) theories as being scalar-tensor theories with a 
gravitational strength coupling to matter, and Palatini theories as being ones with an 
essentially infinite strength coupling to matter (in this sense they also have an infinite 
mass, but the ratio of the mass to the coupling is finite). If one considers a general 
scalar-tensor theory with arbitrary strength coupling (but defined so that the mass 
divided by the coupling is fixed) then one would expect to see a cross-over, for some 
coupling strength, from a behaviour where the macroscopic dynamics are determined, 
to leading order, by the microscopic field equations (and averaging works as one might 
expect), to a behaviour where the average macroscopic dynamics are not described 
by the microscopic field equations. This is precisely the behaviour that was found 
in the study conducted in [T3]. For a relatively weak (e.g. gravitational strength) 
coupling, such as in metric f{R) theories, averaging over small scales works as one 
would generally expect (at least to determine the leading order deviation from GR). 
It should be stressed though that these problems with averaging may re-emerge in 
metric f{R) on larger scales (e.g. those of large scale structure in the Universe). A 
discussion of this scenario is beyond the scope of this article. 

Another point worthy of discussion is that the conclusions of this paper concerning 
the macroscopic behaviour of Palatini /(i?) theories applies generally well to much 
wider class of theories. Specifically, for a given theory on can always write (for a some 
choice of conformal frame) the modified Einstein equations as: 

R". - ^R6''u + A5^, = kT„% + (14) 

where A is some cosmological constant term, Tj^j^ is the energy momentum tensor of 
matter and represents all of the modifications from the standard Einstein equation. 
All of our conclusions of the macroscopic behaviour of Palatini f{R) theories, will then 
apply if, for some choice of A, both terms on the right hand side of Eq. (fT^ vanish 
outside particles. One could then replace kT^^ + by kT^^, and this new energy 
momentum tensor would then only have support where T'^m^ i.e. if T'^nn, describes 
a system of particles separated by vacuum then so does T''^. In any such theory, 
one could, as is done in Ref. [4], deduce the dynamics of particles, and hence the 
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averaged dynamics of a set of particles, simply by considering surface integrals which 
depend only on the form of the metric outside the particles. Outside the particles, 
the field equations in this modified theory reduce to C^i, + AJ'^^ i.e. just vacuum 
GR. In this modified theory, particles in a vacuum would therefore move precisely like 
particles in a vacuum in GR. Since the motion of such particles in the later does not 
dependent on the precise composition of the particles, it is essentially irrelevant that 
we have replaced the original matter energy momentum tensor with a modified one. 
It is straightforward to see that if the motion of particles is equivalent to that in GR 
for one choice of conformal frame, it is equivalent for all choices of conformal frame. 
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